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Abstract 



o: 

^SJ . We present an analysis of the adiabatic approximation to understand when it apphes, in view 



of the recent criticisms and studies for the vahdity of the adiabatic theorem. We point out that 
this approximation is just the leading order of a perturbation series, that holds in a regime of a 
perturbation going to infinity, and so the conditions for its validity can be only obtained going to 



i-Q I higher orders in the expansion and removing secular terms, that is terms that runs to infinity as 



the time increases. In this way, it is always possible to get the exact criteria for the approximation 
to hold. 

PACS numbers: 03.65. Ca, 03.65.Ta, 03.65.Vf 



1marcofrasca@mclink.itl 



I. INTRODUCTION 

The formulation of the adiabatic theorem is as old as quantum mechanics itself [l|. A 
sound mathematical formulation was then provided by Kato 2| with the introduction of 
the idea of a gap condition and finally, Berry |3| showed as phases cannot be disposed away 
so easily. Essentially, it can be stated in a similar way as also applies in mechanics and 
thermodynamics: A system keeps on staying into its initial state provided its evolution 
in time is slow enough. To clearly define the terms of this definition one has to rely on 
the specific framework the system is considered. So, for quantum mechanics, we need to 
specify a Hamiltonian H that evolves in time. Once we are able to solve the corresponding 
eigenvalue problem, we get eigenstates varying in time and if we start from one of such 
eigenstate at the initial time, a slow variation of the Hamiltonian, properly stated through 
a time scale T large enough, maintains the system in this eigenstate properly evolved in 
time and multiplied by some phases. Slowness condition can be properly stated and this 
condition should grant the validity of the given conclusion. 

This idea has been so successful that has found a lot of applications but it appears at the 
foundations of quantum computation |4|-l6| providing a promising way to implement it. So, 
it appeared somehow as a surprise the fact that its condition of validity was not sufficient 
P, |8| . Marzlin and Sanders |7| were able to display an example where the condition for the 
adiabatic theorem to hold just fails opening up the way to a hot debate that is well alive 
yet(9l-l2lj and we can safely say that this matter is not settled yet. 

The aim is to paper is to clarify this question based on a couple of papers appeared in 



the nineties 



22 



23| . The reason of the failure of the condition of validity for the adiabatic 
approximation can be traced back on the very nature of the approximation itself. This 
appears to be just the leading order of a series expansion of the solution of the Schrodinger 
equation. This series, arising as a solution of a differential equation, displays exactly the 
same defects of a more known small perturbation series. These are a typical expansion pa- 
rameter and secular terms. A typical expansion parameter always appears in a perturbation 
series due to the physical nature of the problem at hand. This is fundamental to give a 
mathematical meaning to the series itself. Secular terms represent a major failure in a series 
expansion as they just increase without bounds as the independent variable, time in our 
case, is taken increasingly large. So, before to give a meaning to an expansion as is the 



adiabatic approximation one has to dispose of such terms. Techniques exist to do this but 
we do not discuss them here. Rather, we will show how, going to the next-to-leading and 
next-to-next-to-leading orders, we are able to recover the proper validity conditions for a 
given adiabatic expansion. 

As we will show, in agreement with what is expected in Ref. |18|], secular terms are 
present in a perturbation series independently on resonances being there or not. These kind 
of unbounded contributions can appear also in the more general case as we will show. This 
means that our understanding of the failure of the criteria for the adiabatic approximation 
to hold is quite general. 

The paper is so structured. In secUl] we expose the troubles with the adiabatic ap- 
proximation as stated in literature. In sec JIIII we exploit the deep connections between the 
adiabatic approximation and perturbation theory and how this approximation provides a 
dual expansion to the well-known Dyson series. In sec JIVI we show how one can dispose of 
some examples presented in literature, proving the failure of the standard validity condition 
in the adiabatic approximation, and obtaining the proper validity conditions. Finally, in 
sec|V] conclusions are given. 

II. FAILURE OF THE ADIABATIC APPROXIMATION 



As it is well-known from textbooks [2J], the adiabatic approximation in quantum me- 
chanics shows how a system keeps on staying on the initial state, provided the parameters 
of the Hamiltonian are varied slowly enough. A proper mathematical definition of slow 
evolution should be given for this statement to be considered as a theorem and to have a 
quantitative understanding of the limitations of this kind of approximation. 

The main result can be stated in the following way. Given a quantum system with a 
Hamiltonian H{t) and letting it varying on a finite interval of time t G [0, T], the Schrodinger 



equation admits the following approximate solution for the time evolution operator 22|, |23| 



U{t,0) f^ ^e*^"We-^i'o^*'^"(*'V;i)(«;0| = f/o(t,0) (1) 

n 

provided we are able to solve the eigenvalue problem H{t)\n;t) = En{t)\n,t) and being 
in{t) = {n; t\idt\n; t) the geometrical phase. It is common practice to absorb this latter 
phase into the eigenstates but, as we will show below, this practice should be dismissed. 



This holds if the following condition holds [l|, |2 1 
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\En{t) - Em{t)\ 



{n;t\H{t)\m;t) 



Enii) — Emit) 



<1, 



(2) 



provided a gap exists between the eigenvalues -E„(t) and Em{t). Recently, this condition 
has been proved to be neither sufficient nor necessary marking evidence for the failure of 
the adiabatic approximation even if it seems to hold [7|, |8|. We will discuss several examples 
of this in sec JIVt taken from the literature, and we will also show how to get the right 
adiabaticity condition. 



III. ADIABATIC APPROXIMATION AND PERTURBATION THEORY 



Adiabatic approximation is generally stated for slowly varying Hamiltonians[24]. Typi- 
cally, one has to solve a time-dependent problem like 



H{t)\i;)=thdt 



(3) 



Now, varying t, one has a continuous set of Hamiltonians and if the system was in the 
eigenstate |n; 0) of H{0), the adiabatic theorem states that the system stays in the same 
eigenstate |n; t) of H{t) at later times, except for a phase factor. This can be also seen by 
introducing a time scale T so that our problem to solve can be stated as 



H{t,T)\ij)=thdt 



(4) 



If T is taken large enough, on smaller time scales the system is adiabatic. One can apply a 
change of time variable as r = t/T and write down the above equation as 



TH{T,T)\i;) = ihdr\ 



(5) 



and the solution of this equation for T — )■ oo, at the leading order, is exactly the statement 
of the adiabatic theorem. This can also be said by approximating the unitary evolution of 
the original problem as in eq.([T]). 

The main question about this approximation, arisen in these last years, is its condition 
of validity. Condition ([2]) seems to arise quite naturally from higher order terms of the 
adiabatic approximation. Indeed, given the Hamiltonian 22l. |23| 



H'it) 



Y^ e*[T»W-^'"Wle-s/o'^*'[^™(*')-En(t')](^.t|,;i5^|^.t)|^.0)(n;0|, (6) 



higher order corrections to the unitary evolution of adiabatic approximation are promptly 
obtained and we have the identity 

f/(t,0) = Uo{t,0)Texp (-^ I dt'H'it'yj (7) 

being T the time ordering operator. Eq.(I7]) is equivalent to 

U{t,0) = ^e^^"We-*^o'^*'^40|^.i)(^.0| (8) 

n 
+ 1 \^ e^Me-Tilodt'EUt') f (lt'^ibn{t')-/Ut')]^-j;J^'dt"lEk{t")-E„{t")] ^ 

{k;t'\ihdt>\n;t')\k;t){n;0\ + . . . 

that has exactly the structure of a Dyson series. So, apparently, condition (|5]) should grant 
that the higher order corrections are negligible impeding in this way transitions to other 
states. 



Some years ago, Mostafazadeh 22| proved that the adiabatic series has a parameter 
that characterizes the expansion. This can be easily understood as Hamiltonians represent 
physical systems and so have physical parameters. An expansion for the unitary evolution 
will contain dimensionless ratios of these parameters. Indeed, we proved that the adiabatic 



expansion is dual to small perturbation Dyson series 23|]. This can be seen from eq.(l5]). We 



can write down the following problem for the Schrodinger equation 

XH{t)\^/j)=thdt\^lj) (9) 

being A an arbitrary dimensionless parameter. In the limit of A — ?■ this can be immediately 
solved with a Dyson series 

nt -I ft ft" 

U{t, 0) = / - ^A / dt'H{t') - ^A^ / dt' / dt'dfH{t')H{t") + . . . (10) 

n' Jo n Jq Jq 

or, in more compact notation. 






U{t, 0) = r exp ( --A / dt'H{t') J . (11) 

In the opposite limit A — )■ cxd, one recover eq.(jS]) from eq.([n]) with the dimensional parameter 
T interchanged with the dimensionless parameter A but in the same increasingly large limit. 
So, one has immediately the adiabatic series ([7]) by invoking the proof of the adiabatic the- 
orem. These series are dual each other as the expansion parameter of the Dyson series is 



the inverse of that of the dual series. This appears very clear if we would like to attach a 
meaning to the dual Dyson series in the limit A — > oo. This fact has a serious implication 
for the adiabatic approximation as normally considered in textbooks: Being an ordinary 
perturbation series this is plagued with secularities. A secularity is a term in a perturbation 
expansion that increases without bound as the independent variable is increased. These 
are so dubbed as firstly were recognized in perturbation calculations for celestial mechanics. 
The very existence of such terms in higher order corrections of the adiabatic approximations 
implies immediately that the condition ([2]) cannot grant, as is, the validity of the approxi- 
mation itself. In fact, it is a normal way out from the failure of this condition to take some 
bounds in time for the validity of the adiabatic approximation. This is what was normally 
done when firstly secularities were met in perturbation series but, of course, is not a proper 
solution. The reason is that one has to understand where the problem comes from and to 
solve it correctly. This will be shown explicitly in the next section. 

The main result we will give here is that the adiabatic theorem represents just the leading 
order of a perturbation series. This series displays all the common drawbacks of a standard 
perturbation series making the original condition for the adiabatic theorem to hold, eq.(l2]), 
generally useless for the most common applications. Of course, this could be escaped in 
the original and subsequent formulations of this theorem as a concrete understanding was 



achieved just in recent years 22|,|23||. Removal of these drawbacks recovers a correct condition 



but makes this dependent on the particular problem at hand. 

IV. EXAMPLES 

A. Validating the theoretical analysis 

In order to show how the above stated theorem on dual series works, we consider a well- 
known physical model and compare Taylor expansions of its exact solution with the explicit 
computations obtained from the Dyson series and its dual (the adiabatic expansion). The 
model we are going to consider is the Jaynes- Cummings model, widely used in quantum 



optics and we follow strictly the treatment given in 25(]. It is important to see from this 



example that resonances are not needed to see the appearance of secular terms in the series, 



as correctly advocated in Ref. 18|. The Hamiltonian is 26 1 



Hjc = ua^a + ^(|2)(2| - |1)(1|) + g{\2){lW + |l)(2|a). 



(12) 



This represents a two-level atom coupled with a single mode radiation of frequency u and 
coupling g. 

One can rewrite the Hamiltonian in the interaction picture obtaining 



i7ig = ^(e^^*|2)(l|at + e-^*|l)(2|. 



(13) 



being A = wq — w the detuning that here we take different from for the sake of generality. 
It is immediate to realize that the physical parameter for this model is the ratio ■^. This 
means that any perturbation series and its dual will have this parameter or its inverse as 
a development parameter. The exact solution of the Schrodinger equation in interaction 
picture can be found by looking for a solution in the form 



>/ = ^Ci,„+i(t)|l,r2+ 1) + C2,„(t)|2,n) 



(14) 



being n the photon number. The probability amplitudes are given by 26 1 



Cl,n+l(^) — \ Cl,n+l(0) 



COS 



2 A /fi„t 

. + T^ sin 

2 / fi„ V 2 



2ig^n + 1 . [Vt^t 

-C2,n(0)sm 



i Lr, 



C2,n(i) = { C2,n(0) 



COS 



r2„t 



iA fVtnt 

a:H— 



2ig^/n + 1 



ci,„+i(0)sin 



2 

i Lj-iT' 



3-7A£/2 

(15) 



being f2„ = a/A^ + TZ^ and TZn = 2g\/n + 1 the Rabi frequency. The Dyson series is 



obtained by expanding the above solution in Taylor series of A = -^ giving till second order 



Cl,n+l(t) = { Cl,n+l(0) 



C2,n(t) = { C2,n(0) 



l + ,^(At + z(l-e-^^*)) 



.A 



A 



C2,n(0)(l-e-^^*) + O(A 



(16) 



l-z-(At + z(e^^*-l)) 



4 



77Ci,„+i(0)(e 



•.At 
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0(A=^ 



The secularity in this perturbation series can be easily recognized as this is the term that 
grows without bound in the limit t — )■ oo. In perturbation theory, unless we are not able to 
get rid of the secularity, such a series is not very useful and does not say too much about the 
true solution. Several methods exist to resum secularities as e.g. the renormalization group 
methods described in Ref. 28|. But here, the problem can be easily traced back to the Taylor 



expansion of the functions sin(vT+~?t) in e, having vT+~? = 1 + y + 0{e'^). This can be 
ehminated by simply substituting A with A + ^^ everywhere in the approximate solution 
into the exponentials of eq.f ll6p . The expansion (TT6|) can be immediately obtained through 
the Dyson series ( TTOl) . So, as expected, this series gives an analysis of the Jaynes-Cummings 
model when the detuning A is enough larger than the Rabi frequency 7^„. We have easily 
disposed of the secular term that made the series useless. 

Now, let us repeat the above discussion in the opposite limit with the Rabi frequency 
larger than the detuning. Again, by Taylor expanding the exact solution one has 



Cl,n+l(^) 



Cl,n+l(0) 



COS 



7?.„ \ i , fTZn 



7^„, 

sm I 1 

2 



1 
2A2 



sm 



Tlr 



1 7^„^ . (Un^ 

— -r tsm 1 

2A2 2 V 2 



— tcos 1 

2 V 2 



«'a^ 



-iAt/2 



;i7) 



C2,n(t) = \ C2,n(0) 



7^n 
COS I —t 



X 



sm 



ici,„+i(0) 



sm 



2 J 2A2 



7^. 



sm 



2A2 2 



-tsin 



7^. 



-t 



7^„ 



^t 



Ur 



-tcos 



2 
7^, 



+ 



JAt/2 



2 y 2 V 2 

with the same problem of a secularity at second order. This series can be immediately 
obtained by the adiabatic expansion (dual Dyson series) showing what could seem an unex- 
pected result from the adiabatic approximation, but in complete agreement with the results 
of Ref.|23| and our preceding discussion. To compute the dual Dyson series we need the 
eigenstates and eigenvalues of Hj^. It is easily found that for the eigenvalue gy/n + 1 we 
have the eigenstate 



|a,n,t) = ^(e-*^*|l,n + 1) + |2,n)) 
v2 



and for the eigenvalue —g^n + 1 we have the eigenstate 

1 



|6, n, t) 



v^ 



(|l,n + l) 



jAi 



2,n)). 



(18) 



(19) 



In quantum optics, these are normally dubbed "dressed states" |26|. Geometric phases are 
then easily computed to give 



7a = (a,n,t|z — |a,n,t) = — 



7b = (6,n,t|z — |6,ra,t) = - — . 



(20) 
(21) 



Then, after some algebra using the dressed states computed above, the unitary evolution 
operator is given by. 



cost -^t\ {e-'^*\l,n + l){l,n + 1\ + e'^^\2,n){2,n\) 
«sm| -ytj (e'^^*|l,n+l)(2,n| + e'^*|2,n)(l,n+l|) 



(22) 



that, for \ij{0)) = ci,„+i(0)|l, n + 1) + C2,„(0)|2,n), gives 



\m)i 



+ 



cos ( — ^t j Ci,„+i(0) - i sin ( — ^t ) C2,n(0) 



COS I —^t ] C2,„,(0) - zsin ( — ^t ) ci,„+i(0) 



e"'^*|l,n + 1) 



e*^*|2,n) 



(23) 



2 J ^"'^ ' V 2 

that is the exact form of eqs. (ll7p when higher order terms beyond the leading one are 
neglected, i.e. when A — t- oo, as expected from the results of Ref. J23J and our discussion. 

We improve our computation going to higher orders. Again, using the above expressions 
for the dressed states one gets 

A 



H'it) = --[cos(7^„t)(|l,n + l)(l,r^+l|-|2,r^)(2,r^|) 



(24) 



- ism(n„t){\l,n + l){2,n\ - \2,n){l,n + 1\)] 
so that, the first order correction to the leading order evolution operator Uo{t) of eq. (l22|) is 



given by 



1 



7^. 



Ui{t) = -iUo{t) / dtiH'{ti) =i-sinf — ^t) {e-~'-'\l,n + l){l,n + l\-e'-'\2,n){2,n[p) 
that gives the first order correction 



\5Mt))i = ^\ sin C^t ] (e-t*ci,„+i(0)|l, n + l)- e^t*c2,„(0)|2, n)) 



(26) 



again in agreement with the Taylor expansion as given in eqs. (ll7p . to order j. So, in the 
same way we have at the second order 

U2{t) = -Uo{t) [ dhH'ih) r dhH\t2) = 
Jo Jo 

. 1 



2A2 



sm I 1 1 cos 1 

2 2 V 2 



e^^*|2, n)(l, n + 1| + e"'^*]!, n + l) (2,(^1^) 
- z^t sin (^A (e-'^*|l, n + l){l,n + l\+ e't*|2, n){2, n\) \ 



then, one has 



\52i^{t))i = ^0x2 



1 
2T2 



sml— t -— tcos — t 



(28) 



The agreement with the Taylor expansion as given in eqs.f ll7l) . to order p-, is complete. This 
series is dual to the Dyson series as the development parameter is 1/A. It is clear from these 
computations, by inspecting eq.( l27|) . that the adiabatic expansion produces secular terms 
and these terms impede a correct evaluation of the validity bounds through eq.(|2]). It is also 
clear that imposing bounds on the range of the time variable cannot solve the question but, 
rather, secular terms must be resummed to make the series meaningful and only after this, 
one should establish the range of validity of the adiabatic approximation. 

B. More common examples 

In order to give one of these examples of failure, we consider the one discussed in |8|, Il8| . 
This example has the property to be exactly solvable and the adiabatic approximation can be 
immediately compared with the exact solution. The Hamiltonian, firstly due to Schwinger 



27|,is 

H = (o"a;Sin6'cos(X't + cTy sin 6' sin wt + a^ cos 9). (29) 

Instantaneous eigenvectors and eigenvalues are immediately obtained as 

huo , , / e~*^*cos| \ , , / e"*^*sin| , ^ ^ 

Eo,i = T^, |0;t)=( ; , |l;t)= I I, (30) 

e2'^sm| / \ — e 2 cos I 

and the geometrical phases are given by 70,1 = ±^cos^. So, the adiabatic approximation 
provides the unitary evolution 

that is 



. ( cos ^ut + i cos 9 sin ^ut) e ^'a isin6'sin^e ^'2 

Uo{t, 0) = P 2 ' ^ 2 (32) 

t sm 9 sm ^e 2 (cos ^ut — i cos 9 sm ^wtj e 2 



10 



being u = uq + oo cos 6. The exact solution takes the form 



a 



(cos ^ + i '^+^^^^'^g sin ^)e-^f i^^^a^M gin ^e-^1^ 

?7(t,0)= ^2 " _ /^ _ '^ ' _ , , (33) 

^!£o|ine sin ^e^^ (cos ^ - i '^+'^iL^"'^^ sin ^)e^^ 



being w = a/wq + i^^ + 2a;ou; cos Q. Now, we can recover eg. (1521) from the exact solution 
quite easily with the condition wq ^ w. In this case, cJ ~ cD and it is consistent with the 
condition for adiabatic approximation that is easily obtained from ([2]) to be cuq ^ i^ sin Q. 
So, we realize that the adiabatic approximation is just a standard series expansion as any 
other perturbation expansion but we discuss this below. Here, we just note that adiabaticity 
condition proves to be not sufficient to grant the applicability of this approximation. This 
can be immediately realized if we try to apply it on the exact unitary evolution to get 
higher order corrections. These cannot be neglected if we take the extreme of the interval of 
time evolution of the adiabatic approximation T large enough and so our condition fails to 
grant applicability. Now, using eigenvalues, eigenvectors and the geometric phases already 
computed, we are able to get the Hamiltonian for higher order corrections 

i/' = -l;ia;sinee-^("'°+'^^°^^)*|0;0)(l;0| + c.c. (34) 

We see that when the resonance condition wq + i^cos^ f^ is met, a secularity is immediately 
produced causing the failure of the standard adiabaticity condition. Due to the presence 
of these terms, the adiabaticity condition cannot be neither sufficient nor necessary unless 
these terms are properly summed away. Particularly, for this example, even if the resonance 
condition is not met, going to higher orders will give rise to powers of time variable making 
the adiabatic series more and more unreliable. At resonance one gets 

H' = sin^ 9a z H sin 9 cos 9ax (35) 

and the secular contribution is then 

- i / dt'H'{t') = i — sin ^(sin 9az - cos ^a^)t. (36) 

Jo 2 

The existence of such a term makes useless the adiabaticity condition. But when we are off 
the resonance, the situation is the same as the condition for adiabatic evolution just fails. 
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In fact, the second order term for the adiabatic evolution can be easily evaluated as 



dt'H\t') / dt"H'{t 



1 

~4 

+ e 
1 



hWsm^e [ dt' 

Jo JO 

i{u)(,+u}cose){t'-t") 



dt" 



-i(^o+-cose)(i'-t")|o.o)(0;(|^7) 



1;0)(1;0 
t 



^nVsin^^- 

4 ujq + oj cos 



[|0;0)(0;0|-|l;0)(l;0|] + n.s.t 



with n.s.t. meaning "not secular terms". So, for this kind of model, the next-next-to-leading 
order produces a secularity in any case that makes the adiabatic series fail. This term must 
be resummed to get a fine adiabatic expansion. We just note that we can read off the 
adiabatic condition in the form \uj sin 9 / [uq + u cos 9)\ ^ \usm.9 /uj\ ^ 1, in agreement with 
the conclusions given in 18|]. This is also the expansion parameter of the series. When t is 
increased enough, the series fails and the adiabaticity condition becomes useless. 

Another simple and largely known example is the one given in 19|. The Hamiltonian is 



H = Cr^ — V COSUJotcTx- 



(38) 



Here and in the following ax, Cy, az are the Pauli matrices. This model is widely studied 
and it is known that it undergoes suppression of tunneling at the leading order but the 
corresponding condition is not exact being lost by higher order corrections. All this is 
computed through the adiabatic series exposed in sec JIIII 29l . |30| but see also 3l| for a 
different approach. The essence of the approach to recover the adiabatic approximation, as 
already stated above, is the elimination of the secularities into the adiabatic perturbation 
series. In ref.[19| it is claimed that rotating wave approximation should apply but this is 
not consistent with the adiabatic approximation, as we are going to show in a moment. In 
fact, it is not convenient to apply directly the adiabatic approximation to this Hamiltonian. 
Rather, we do a change of representation going to the interaction picture. This will give 



Hf = — V cos Uota^e 



iiuzt 



and so we get the eigenstates and eigenvalues 



£"0,1 = ±V COSUot, 



|0;^)^=^ 



v^ - 



eh 



et 



i^^*>'=7! 



e-^* 



(39) 



(40) 



and for the geometric phases we get 70,1 = =Ff • So, in the interaction picture, the unitary 
adiabatic evolution is given by 



[/o/(t,0) 



e 2h^^e ^"0 |0;t)7 /(0;0| + e2fi"e'"^o "| 



l;t)//(l;0| 



(41) 
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and finally 

?7o(t,0) = e^^'^^U/(t,0) (42) 

that gives 

f/o(t,0)=e"^^""""*|0;t)(0;0| + e^^"""''*|l;t)(l;0| (43) 

being 

1 / 1 \ 1 / 1 \ 

(44) 





that are independent on time and are eigenstates of ax- In this way we can immediately 
compute the Hamiltonian for higher corrections obtaining 

H = — (T^e'*"o "=. (45) 

With this Hamiltonian we can evaluate the correction to be neglected for the adiabatic 
approximation to hold. Indeed, using the identity 

^iJ_sina;ot<x. = J^ Jn ( ^^ 6^""°*"^ (46) 

being J„ the Bessel functions of integer order, we recognize that 

and this constant term will give rise to a secularity producing a failure into the adiabaticity 
condition. But this term can be just resummed away to give a proper unitary evolution 



29 



3Q| 



U{t,0) ^ e^^"""°*"^e^^^''(^)*"^ (48) 

and the approximation holds provided eJn{V/huio)/fi^o ^ 1 for any n > 0. This is again a 
strong coupling expansion as this condition is granted, for any n, provided V ^ JtjjJq also at 
resonance and this rules out rotating wave approximation. But we turn our attention to the 
adiabatic theorem, we now realize immediately that, after resummation of the secularities, 
this holds unconditionally while the usual adiabaticity condition, Vloo\ sinwotj/e^ ^ 1 (see 
[l9|), just fails. So, to verify how and when the adiabatic theorem holds we need to go to 
the next-to-leading order, at least, and resum secular terms. 
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V. CONCLUSIONS 

We have shown the origin of the difficulties with the standard vahdity condition of the 
adiabatic approximation. This just arises from the perturbative nature of this approximation 
that makes enter into play some singular behaviors of the expansion. Disposing of these 
terms puts out the right conditions for the validity of the adiabatic theorem. Besides, a 
clear understanding of the expansion parameter of the series should be achieved to state 
properly a validity condition. 

As a side result, we are able to produce a perturbation series for a strongly perturbed 
Schrodinger equation. This is a quite general approach to manage equations with a term 
going formally to infinity. 

Our hope is that this should give rise to a new understanding of this old approximation 
providing some new insights in fields of physics where applications are at hand. 
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